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Abstract. Various procedures for expressing the multipolar expansion of the 
electromagnetic field are considered with application to the calculation of the radiated 
power. Some results from the literature are discussed and perspective of developing 
the subject is pointed out. 



1. Introduction 

The multipole expansion of the electromagnetic field is an usefull tool in physics and 
can be found in any book on electrodynamics or on the theory of atomic and nuclear 
transitions. In most of the textbooks (see, e.g., [UEj) a full and systematic treatment is 
given in spherical coordinates, while for Cartesian coordinates the problem is presented 
fully only for the static case; the dynamic case is given only for the lowest order 
multipoles. A general procedure for the reduction of the multipole tensors represented 
by Cartesian coordinate components to fully symmetric traceless ones in the static case 
is given in [HI Ej and it is generalized in to the dynamic case. This method is applied 
in Hj to the radiation field. 

In the present paper different procedures are applied for calculating the total power 
radiated by a confined system of charges. A first one represents the traditional method 
of expressing the field expansions by the multipole Cartesian tensors, applying finally 
the reduction of these tensors. The second one uses the reduction technique done in 
but more systematically considered here. The results of jHllZj are analysed and compared 
to those of other methods from the literature [HI ED]. The advantages of the Cartesian 
coordinates are emphasized firstly by the simplicity of the formalism : only algebraic 
manipulations and combinatorics are implied, no special functions being required. 
Secondly, the procedure initiated in (Hj and used here leads to a nontrivial grouping 
of different multipolar contributions standing out the toroidal multipole contributions 
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In section 2 some basic formulae for multipole expansions are presented. Section 3 
deals with the radiation field as well as with the expression and expansion of the total 
radiated power. In Section 4 the procedure of reduction of a tensor to a symmetric 
traceless one is given. The total radiated power is then treated in Section 5 using these 
reduced moments and a comparison with literature is made. Section 6 presents, in a 
more systematic and concise way, the results from [SJ U\ ■ Then the total radiated power 
is expressed in Section 7 by the transformed moments. The conclusions are given in 
Section 8. In Appendix the proof of some formulas used throughout this paper is given, 
as well as the reduction scheme and a justification of the reduction procedure using the 
charge and current expansions. 



2. Basic Formulae for the Multipole Expansions 



Let us consider charge p(r,t) and current j(r,t) distributions having supports included 
in a finite domain T>. Choosing the origin O of the Cartesian coordinates in V, and 
using the notation d for the orthogonal unit vectors along the axes, the retarded scalar 
and vector potentials at a point outside V, r 



$(r,t) 



p(£,t-R/c) 



Po fj(£,t-R/c) 



d 3 e 



4tt£ J R 4vr J R 

where R = r — £. The Taylor series expansion of the function f(R) is 
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and a n is the 77-fold tensorial product a <g> . . 

Denoting by T n a nth order tensor, A ( - n - > ||B < - m - ) is a \n — m|th order tensor with the 
components: 



(A^HB^. 



in-mjl— jm-Bji---j m , n > m 

Bn-,...rj , n — m 

31 Jn Jl Jn 
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By applying the formula for the Taylor series expansion to the scalar potential we get: 
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pM being the n-th order electric multipole tensor. 
For the vector potential we obtain the expression: 



, pW(t) = / r"p(r,t)d 3 x, 
v 



(2) 



**h n! 



^+ 1 ){t-r/c) 



^E — r e A"-^ 



Ph-i n i(ta) 



t = t — r/c 



(3) 



Some Aspects of the Electromagnetic Multipole Expansions 



In the previous equation the magnetic multipole tensor was introduced by its Cartesian 
components: 

r> 



3. The Radiation Field 



The formula for the power radiated by a charged system described by the charge p and 
current j densities with supports included in a finite domain T> is well known [2]: 



r 

Hoc 



d . , . 
v x —A rad (r,t) 



(4) 



Here, u = r/r and dP/df2 is related to the flow of the energy detected in the observation 
point r at large distance r compared with the dimensions of the given charged system. 
The vector A ra d is obtained from the retarded potential ((J) by retaining only the 
dominant terms at large distances. 

In the following the expansion of A ra d is derived [UIE]. Starting from: 

A(r, t) « f ± £ Li VI ^(f - r/c) 
4vr r ^ n! 

« j^ei- jr t-^-di, ■ ■ ■ d in n iv .. ini (t - r/c) 
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and considering : 

9n ■■■d i J(t-r/c) 



-^^■■■^f(t-r/c) + 0(l/r) 



one obtains the following expression for the part of A contributing to the radiation: 
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with the notation: 



(n+l) = U (n+1) 

So, considering the angular distribution of the radiation given by equation (JH) and 
applying the expansion of A rad , one gets: 
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Finally the following result for the angular distribution of the radiation is obtained: 
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d m+1 
X dt m+1 ^ jl '" jm ' k ' 

n>0,m>0 ' t -" i -<- 

In order to calculate the total radiated power the formula J = 4tt < J(y) > is 
used, where < f{u) >= (l/47r) / f{u)dVL{u). The result is: 



j=—{ r 1 

^7rfc C- n >0,m>0 "'■'"'■ C 



d n+l d m+l 
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In the previous equation it is usefull an averaging formula introduced in jHJ d and 
justified in jH]: 

When working with the expansion (JJJ) of the radiated power, one must have a strict 
criterion regarding the comparison of the different terms contributions. This criterion 
can be easily obtained if one refers, particularly, to monochromatic sources. If the 
possibility to represent any type of variation in time as a superposition of functions 
corresponding to the monochromatic sources is taken into account, the conclusions will 
be general. In the case of the variation characterized by the pulsation u, the term from 
the expansion indexed with the pair (n,m) contributes with an order of magnitude equal 
to (d/\) n+rn to the radiated power. Since the expression of the radiated power comes 
from the values of the field in the zone where A >> d, it is obvious that the order of 
magnitude of a term is given by the sum n + m. Therefore, a consistent way of using 
a finite number of terms in the expansion Ffy, terms that accurately characterize the 
radiation up to a given order M, is to retain all the terms corresponding to the values 
n + m between and M. Apparently, in literature such a procedure was not consistently 
used everywhere, as it will be shown later on. 

So, considering j( M ^> the M order term from the expansion (JJJ), we have 

J-(M) _ j-(n.m) 
n+m=M 

4ne c 3 J (n ' m) = 

<HiA + a - (""W K>SS?)>- (») 

For particular cases, 

j(0) _ J(P,0) _ j(0,l) + j(l,0) _ q 

J m = j(o,2) + j% ) + J (i,i) i t7 (3) = 0j 
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where the fact that the averaged terms with n + m odd are zero is considered. 
Below we will present the first terms of the expansion: 
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For j( 2,2 \ JT^ 1 ' 3 ', J'^°' 4: \ the results of the contractions with the 5— tensors will be given 
bellow by a simpler method. 

The expressions for J7"'- ' '', J^ ,2 \ J^ 2 ' ** are given by Bellotti and Bornatici 

in They go further, introducing the reduced multipole moments and finding a new 
term, as we will see in the following. 

Sometimes it is easier to use the magnetic moments defined as: 
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It can be shown that instead of the expansion (JHJ) one can use an expansion obtained 
from this one by performing the substitution: 



t i i\...i n y £i n i n _ 1 feMi 1 ...j n _ 2 fc + P«i...i n 
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The result for the vector potential is given in [H |Hj : 
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The transformation (fT2"j) is used in the equations (|1U|): 
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Here the fact that the contractions of symmetric and antisymmetric pairs of indices 
cancel is considered. Because < UiUj >= (1/ 3)<5^-, < v ix Vi 2 > ea^rhk Pi 2 i= 
(l/3)eu lk rh k P hi = etc and Ea^e^i = 28 kh so what is left is: 
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Here, again, there is a notation 

|\K 3,1) — £■ M 
ij tp s > tps- 

4ne c 3 J^ = 4ne Q c 3 J^ 
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The terms where the symmetric-antisymmetric tensor contractions are present cancel 
again, giving the final result for J^ ' 4 **: 



4n£ c 3 J^ = Ane c 3 J^ 



24c 4 



4. Symmetrising and "Detracing" the Tensors 

In this section the symmetrisation and detracing method specific for the electromagnetic 
moments is presented. The procedure and notations from |lj and [S| are used, as well as 
a theorem introduced by Applequist in [Tlj. The first step is to symmetrise the magnetic 
tensors M^, already symmetric in the first n — 1 indices: 



M {sym)h ... in = -[M h .. An + M in ... h + M,., , . + . . . M tl . 



•Ini^n— 1 J 
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where represents the sum over the independent terms only, from all the permutations 
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of the n indices. It can also be written as: 
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where is the component which does not have the i\ index, while: 

N^ n,1) =£■ M- ■ 
is a n — 1 order tensor. Further, the correspondence: 

T (n) -> M[T {n) ] 
is introduced, where TV^T^] is a n — 1 order tensor: 

•A^T 1 - ^]i!...i„_i — £i„-ips~T i 1 ...i n _ 2ps 

and A/" fc [M ( - n - ) ] = U^ n ' k ' is a tensor of rank n — k. Here are some examples: 
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The following relations are also to be considered: 

7V 2fc [M (n) ] = N (n;2fc) = fclL^ /" ^2^n-2fc x jtfg = (_l)fc|V| (n;fc) , 

n + D 

^2fc+lr M (n)j = N (n;2fc+1) = (-1)*^ f £2fc£n-2fc-l x x j)d 3 £ = (-l) fc |\| ((n ' 1);fc) , 

fc = 0,1,... 

where j^ n ---^ k i s the tensor obtained from the contraction of k pairs of indices. 
Particularly, the tensors introduced by equations (fT3j) and f)15jl are 

N i = N?' 1) J N^-^Ng^. 

For the reduction of a totally symmetric tensor to a traceless one we write as in 

£>ii...i n — 5'u...i n — J]] 3iii2-h-[S^ ]i3,i4...i„ 
D(i) 

where S is a totally symmetric tensor and AfS 1 ^] is a totally symmetric tensor of rank 
n — 2. Further the following notations are used: 



A[M (n) ] = A (n ~ 2) , A[P (n) ] = n (n " 



-2) 



Moreover, sometimes, for writing some equations in a simpler form, we will also use the 
notation A[T^] for an arbitrary tensor but by this notation we will suppose that 
the symmetrization of ~P n ' is implied i.e. 

A[tW]=A[t£J- 

In ^3] a general procedure of detracing a symmetric tensor is presented. In our 
case the result of this procedure may be written as: 

A[ fl w 1> . A =E ( - 1) 7 a, :,:- 2m) " x 

[2n — ly.lm 



m=l 

13 14 ■ ■ ■ w l2m-l I2m'- J i2m+1 
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Using the definitions and notations presented above it is easy to obtain the following 
useful results: 

• detracing the electric moment tensors: 

II = A[P( 2 )]: n = ^P« (16) 

n (i) = A[P (3)] : Ui = ^P qqi (17) 
n( 2 ) = A[PW] : % = ±P qqij - ±f qql A, (18) 
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n (3) = A[ p(5)] . n ,, fc = lp 9 ^ fc - TI 1_ £ ««p«b* 
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• detracing the magnetic moments tensors: 
A = A[Mgj = 

A (1) = A[Mg)j : A, = ±M qqi 
A (2) = A[M«] : A, J = 4(M WJ + M, 



2g \ "WJ ' 

A (3) = A[M jgj : A iifc =i^M 93iiA 1 
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5. The Total Radiated Power in Terms of Reduced Moments 

The usual procedure is to emphasize, in the expression of J~, the reduced tensors by 
utilizing the reduction relations given in the previous section. The well-known notations 
for dipole moments are used: 

Pi =Pi'P, Mi = mi : rh 

The static expressions of the reduced tensors of electric and magnetic polarizations are 
V {n) and M {n) with : 
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Below there are the results of the reduction of the first terms in the expansion of J~, 
which correspond to the limitation to the order (<i/A) 4 : 
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3c z 4 o 
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The toroid dipole moment fTJ was introduced, with the definition: 

4 o 

There was also used the simplified notation Nj = N 4 - 2,1 ' ) , while IT^ 2 ) is a traceless tensor 
(that is "detraced"): 

111 

The sum of terms from the equation (J24|) represents the total radiated power expanded 
until the fourth order with respect to d/X, that is: 

n+m<4 

It is verified the relation: 

2-2 4 - 1 1 - .:. 1 .:. .:. 2 /■■ 1 -\ 2 

~ 3 P ~ 3^< T, +— A N,N, -— N, n, + s n, ni = 3 (p - ? T) 

By eliminating the higher order derivatives, we may write: 

1 3- 1 ~ 3 ~ 3/1- 1 ~ \ 

— ^ikPik — —z^ikPik = —^ikPik — —^-ikPik = — -Nj& — — Ilifc Pifc 
15 40 15 40 10 \9 4 / 

(3 1) nr^n _ ^ _____ ____ 

with the notation: Nik = N) fc ' . The so-called quadrupole toroid moment [TTJ [T21 UHl US] 
is also defined by: 

Tjk = — Nit — — rh/o 

Explicitly, the two toroid moments are written as: 

v 

= \j '[{*• m - etm - ^ / 

V V 



but since: 



/ z%f>a 3 z = -J ■ jd 3 z = jj. (£%)a 3 z = J [2(1 ■ m + 

V V V V 



it results: 
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In a similar way: 
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After regrouping the terms and considering the above presented relations, one gets: 

, ,^ 2 /-••. 1 I-v 2 i 1 1 ji ... 

^J^ = -(p--l) + -m + — VvPv - mTljVli 



4 1 A" (4,1) 4 ::; \ 1 2 •■ •■ 

+ ^ m ^ M ^ - w^Pi [ H m ~ 5 n ^ J + MijMij +-^V ijk V ijk (25) 

If the expansion of the radiated power up to the second order with respect to d/X is 
considered, only the following sum will be kept: 

J(2) = j(0,0) + + 2J (0,2) 

■■l-> 2 --2 4 - 1 

3 P + 3^ m ~ 3^ T * +20^ ^ 



47T£ C 3 

This is, actually, a result obtained by Belloti and Bornatici [Oj. They observe that, 
compared to the expression for the dipolar magnetic- quadrupolar electric radiation from 
the books of Jackson and Landau pH |2] as well as from many other electrodynamics 
books (including the one by C. Vrejoiu @]), this result contains a supplementary term 
represented by the contribution of the vector T. We point out that in QjJ |2] the goal is 
only to calculate the isolated contributions of the electric dipole, electric 4-pole and 
magnetic dipole to the total radiated power without regarding it as a result of an 
expansion. If one associates a multipole moment to an elementary system, one obtains 
an isolated contribution of this multipole, but when one considers a composite system, 
all multipoles, giving contributions of the same order of magnitude, must be considered. 
This is the reason for which the result from PJ Ej should not to be considered erroneous. 

Without knowing, probably, the results published ever since 1974 by Dubovik et 
al. JI] and all the other publications refering to toroid moments, they consider this 
term as being something new. On the other hand, Bellotti and Bornatici give a correct 
result, if one considers how far they go with the expansion of the radiated power. 

The impression is that in general in literature similar results contain a series of 
confusions due to the inconsistent application of the criteria according to which different 
terms of this expansion must be compared. For instance, in 13J a result which is 
supposed to represent the expansion including the fourth order terms is presented: 

2 /•• 1 '^\ 2 2-2 
^Tie c 3 J Dubovik = - f p - — J J + —m 

+ 2~o^2 ^ ik ) {^ ik 

+ {ViikP ijk + k jk k jk ) +■■■ (26) 

Appart from the mistake (of course, typing error) that is derived only three times, 
there are certain things that must be observed: 

• the term Aiij k A4ij k from this expression is of the sixth order (with the 1/A 
criterion); 

• the term T ik T ik is also of the sixth order; 
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• these terms should not be present if in the expansion of J the terms up to the sixth 
order are not consiered; 

• the term M.^M.^ is missing, as well as other terms that should be present, at least 
up to the fourth order; 

• the sign three points " . . . " at the end of the expression is losing its usual sense 
when a finite number of terms is considered in a series expansion. 

In some more recent papers, as for example in [TU] . it is claimed that, as an 
application of more general formulae, the expansion up to the order 1/c 5 is given (that is, 
exactly how it is considered in 9J). The 1/c criterion is correlated with the wavelength 
criterion if one takes into account the powers of the parameter 1/c together with the 
orders of the partial temporal derivatives in the expansion of the radiated power. The 
expression given in JU] (written with our notations) is: 

47re c 3 J R „ v = 1^-1 fj + J^m k M qqk + fijfij MijMij (27) 

where, in order to identify it with the equation (4.12) from [TO], one must consider the 
definitions of the quadrupol moments and the identity: 

this last parameter being the one used in the equation (|?7|). 

It is easy to observe that if one tries to interpret equation (|2*7jl as the expansion 
of J up to the second order, the fourth order term T 2 is present without justification. 
On the other hand, if the same equation is considered as an expansion up to the fourth 
order, many terms are missing. 



6. Reduction of the Multipole Tensors and Gauge Invariance 



In this section the results from U\ will be presented in a more systematic and 
concise way. A different explanation of these results is given in Appendix A using 
the formalism of Dirac's 5— function. These results present the possibility of expressing 
the electromagnetic potentials, as well as the field, exclusively in terms of reduced 
moments, that is moments represented by reduced tensors (symmetric and traceless). 
This procedure allows, particularly, to express the radiated power in a very simple 
general form, as it will be shown below. 

Let us consider the expansion of the potentials with the moments (see equation 
(HQ)) and P (ra) (see equation ©): 



A(r,t) 



-1 



i n-1 



in— 1 1 



Aix 



£ 

71 = 1 



nl 
D n -\ 



rn-ll 



nl 



-M (n) (t-r/c 
r/c) 



ipW (( 
r 



(28) 
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-p(«)(t- r / c )- 



47T50 n\ 



(29) 



The basic idea is that the tensors and can be replaced by reduced tensors. 
This leads to a modification of some of the inferior order moments, such that finally 
all the moments up to (and including) the nth order moment must be expressed by 
reduced tensors. Not all of them, though, reduce to the static expressions M. and V . 
The final results for these moments are M and P. Formally, as a result of this procedure, 
the expressions for the potentials expansions up to the nth order will be given by the 
equations (J2EI) and ((201) witn the substitutions P — > P and M — > M. 

As a basis for this procedure one could take the following observations resuming 
systematically the results from [H]: 

• Let be a -type tensor i.e. 

L is fully symmetric in the first n — 1 indices, Lj 1 .„j n _ 1 i k = if k < n — 1. (30) 
The symmetrisation of L is realised with the relation: 

1 71—1 

\-(sym)ii...i n = Ui...i n — ~~ E £ i\i n q-^'i 1 ...i n ^ 1 q[^ ^] (31) 

71 A=l 

with: 

hf. . [I = £ . I 

In this case: 

a) the substitution: 

i n— 1 

M n ..,„ — > M (i)il .. in = M,,. in - - £ e^^X'k-xJ 1 -^] ( 32 ) 

n A=l 

produces changes of the potentials which, up to a gauge transformation, are 
compensated by the following transformation: 

p(n-i) _^ p(n-i) _ Ii_IXf[L {n) ]; (j) (33) 

b) the substitution: 

i n— 1 

P«.., n — P(L)n.., n - - E ^^l^J^] (34) 

H A=l 

produces changes of the potentials which are compensated by the following 
transformation: 

M("-!) — ► M (n " 1} + ^-^Af[L (n) ]; (II) (35) 

• Let be a tensor and the detracing operation for it: 

^>i±...i n — $ii...i n ~ E/ ^Wi^-iz—in ^] (36) 

Then: 
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c) the substitution: 

M n ..,„ — > M^j, - £ 5 nJ2 A l3 ..,JS (n) ] (37) 

D(i) 

produces changes of the potentials which, up to a gauge transformation, are 
compensated by the following transformation: 

M( n - 2 ) — ► M (n_2) + ^-^A[S (n) ]; (///) (38) 

d) the substitution: 

Pn...i n — P { S)n.., n ~ E <WW~*» [S W ] (39) 

produces changes of the potentials which, up to a gauge transformation, are 
compensated by the following transformation: 

p(n-2) _^ p(n-2) + lZ_2^r S (»)]_ tjy} ( 4Q ) 

These four transformation relations of the electromagnetic potentials are sufficient for 
the development of a scheme in which the replacement of the multipole moments tensors 
by symmetric and traceless tensors is presented. 

Such a scheme, valid for expansion up to the fourth order with respect to d/X, can 
be found in Appendix B with results written in the Tables from Appendix C. It can 
be easily continued for higher orders. It is possible also to give results for the reduced 

~(n) — (n) 

tensors P and M for arbitrary n which will be given elsewhere ^B] • 



7. Radiated Power Expressed by Transformed Moments 



By applying the transformations (jl2j) to the expansion of the vector potential of the 
radiated field, one gets: 

.. 1 OO 1 

A rad (r,t) = !±-^—[vM\\M%\t-r/c)] x „ 
47r r ~y n\c n 

/ioclg 1 (n _ 1} ( } 
4vr r ^ n\c nl 11 ' n y 1 n 
The radiated power (angular distribution) will be given by: 

1 



OO OO 

EE 

n=l m=l 



\ rn \ r n+m 



rv.mic 



^IIM^) (u^\\M^ 



+ EE 

n=l m=l 



n\m\c n+m 



u \\",n+l ll r ,m+l 



oo „ 

+ ee s -^{K- 1 »iim!;i, 



n=l m=l 
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Considering the procedure of reduction of the moments tensors from the expansion of 
the vector potential applied up to the /ith order for the magnetic and to the eth order 
for the electric moments, the sum of the terms from the expansion of the radiated power 
which contain exclusively magnetic moments reduced up to \i and electric moments up 
to e is: 



oo oo 



EE 

n=l m- 



n\m\c n+m 



(n) 
,n+l 



(m) 
m+1 



- (>)||M (n5 

oo oo 

+ EE 



,n+l 

c 2 



u {m) \\M 



(in) 
,m+l 



n=l m=l 



\ rn \ r n+m 



filmic 



,>-i)|ip (n) 

U ll r ,n+l 



i/ (m_1) ||P 



(m) 
,m+l 



oo oo 

+ EE 



n=l m- 



\ n\m\c n+m 



(n) 
,n+l 



(n) 
,n+l 



1/ X [i/^-^IIP 



(m) 
m+1 



(41) 



It is easy to understand that the above sum cannot be identified with the expansion 
of the radiated power containing the magnetic moments M^, k < fi and the electric 
moments P^, I < e. This happens because some of the magnetic moments "* contain 
tensorial expressions built with magnetic and electric tensors of superior orders in k and 
similar for the electric case. But, once J^ £ is settled, the superior orders contributions 
can be eliminated in order to obtain the correct expansion with the d/X criterion. 

Returning to the expression of the total radiated power as a function of the magnetic 
moments and considering the expansion Q, it results that this one corresponds to 
the expansion (|4*T]) for fi = M and e = M + 1 when only the contributions of order not 
higher than M are retained. 

According to the results from :7], expressing the total radiated power is easier if 
one uses the averaging formula (jH)) and the symmetric and traceless character of the 
reduced tensors. Hence, the following properties are valid: 

a) let two symmetric traceless tensors A*™) and B*™) and their averaged contraction: 

( (u k \\A^) || (u k '\\B^) )=<v n ... v ih v h . . . ^A l ,.. ifcJfc+1 ... Jn B J1 ..,,, fc , +i .., m > 

This average is different from zero only for products 5i j k with p 
q = 1, . . . , k' and the following relation is valid: 

k\ 



k and 



((^||AW)||(^'||B( m ))) 



^(n)i ig(m) 



k'k 



{2k+l)\\ 

b) The terms of the last sums, containing mixed moments products, in equation 
(JUJ), give contributions of the type: 



< z/j x . . . Vi n _ x Vj x ■ ■ ■ Vj m -\Vp > c"j nP<? Aj 1 j n Bj^ 



■Jm-iq 
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but all the terms from the 5 products, representing the averages of the v components 
products, contain either 5 ikP or S P j p k = 1, . . . , n — 1, / = 1, . . . , m — 1 such that, because 
of £i nPq and of the symmetry of A and B, the result is zero. Using these properties one 
can write the result: 

>-J [IE 



4ne c 3 



E 



n + 1 



\ nn\(2n + l)\\c 



2n 



M 



(n) 



M 



(n) 



.n+l\\' ,n+l 



+ E 



n 



,n+l 



^nn!(2n+ l)!!c 2 ™- 2 V ,n+l1 

Let us consider now different expressions of J^ M \ calculated this time starting from 
the formula given by the equation written above. As settled, for a given value of M one 
must consider this formula for /i — M, e — M + 1. 

/j, — 0, e — 1, only the electric dipole moment has a contribution: 

JW : (^e) = (2,3); in this case the results from the table in Appendix A are used, 
considering in the end only the terms corresponding to the order M: 



M„M„ .f-PiPi + _ PyPy 



where the index of the bracket is the maximum order in d/X which has to be retained 
in this bracket. It follows that: 



47T£ C 3 J r(2) 

JW: (/i,e) = (4,5) 



2-3 2-^2 4 £ 1 ^ ^ 
3 P + 3? m ~ 3^' + 20^« 



— (MH -A -N )i 

3c 2 V 6c 2 18c 2 h 



+ 
+ 



1 ... 1 -i 1 £< 4 - 2 > 
26? ( ^ + 4? A "24? N 



... l 7- 1 £< 4 - 2 > 
^ + 4? A -24? N 



2 ••■ ••■ 1 

9A5^ MijkMijk + 18144? 



MijkiMijki 



+ I(p 



-— ^— X[N (4,1) 1 + — X [n (3) i + — 
32c 4 y 20c 4 J + 96c 4 



.7.(4,3) 

N 



(P- 3 T 



1 SJ/mJ 20c 4 



+ 



+ 



32c 4 



20c 2 



1 



;VijkiVijki + 



945c 4 
1 



1 - (4,3) 

a rn^i + _ -n )i 

[ J 96c 4 ; 



c 
V 



ijklqP 'ijklq r 
J 4 



2 5 x 3 4 x 7 J J 4 x 3 3 x 5 3 x 77 j 4 

Only the fourth order terms must be kept from this expansion. This could be done by 
eliminating the useless terms, but, for obvious reasons, a detailed description is prefered 
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^e^J^ = "-('p-- 2 J 



here. The table with reduced moments is examined and the terms corresponding to the 
considered approximation are retained from the expressions representing the squares of 
these tensors. The result is: 

2 

3 V" c 2 

A ( 1 •■ 1 1 •■(4,3)\ 

3 F \ 32c 4 [ symi 20c 4 1 J 96c 4 J. 

1 _L_ ... Jl 2 ^ « 

+ 20? ~10? ^ * J + 945? ijk ijk 

2-^2 4 ( 1 ■■ 1 ::(3,2)\ 1 

+ 3? m + 3? m? 6? A -18? N , + 20?- M -- M - 



Because 



N (4 ' 1} 1 

qqi 



and 



v 



the last expression of is the same as the expression given in equation (|23j) . 



8. Conclusions 



The multipole expansions in Cartesian coordinates are not largely treated in literature 
or, if they are, many inaccuracies and uncertainties appear. In papers like [5] and [Zj 
the expansion in Cartesian coordinates of the electromagnetic field is presented and 
applied to radiation. Their results were presented here as well, in a more systematic 
and concise way. They were compared to other results from literature, usually obtained 
by expansion in spherical coordinates. It was underlined the fact that the problem is 
not accurately treated everywhere, this fact leading to some confusions. 

The main difficulty in the case of Cartesian coordinates is the procedure of reduction 
of the nth-order multipole tensors to symmetric traceless ones to obtain in this way a 
description of multipoles in terms of irreducible rotation group representations. This 
procedure was presented here, in Chapter 3. The transformations implied in the 
reduction procedure were defined such that the electromagnetic potentials are altered 
only by gauge transformations. This implied a specific feature of the dynamic case: 
the redefinitions of the multipole tensors in the lower n < N orders, induced by the 
reduction of tensors in a given order N. 

The following problems are left for a future discussion and they could be some 
possible research topics: 

• completely systemizing all types of moments grouped, in the different orders, in 
reduced moments; 

• defining the singular distributions associated to toroid moments and their physical 
meaning; 
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• study of the implication of the existence of interactions associated to toroid 
moments and the setting, according to correct criteria, of the different terms 
contributions; 

• making symbolic programs which realize a reduction scheme of the type presented 
in Appendix for the general case. 



Appendix A. Equivalent multipole expansions 



Following the basic ideas from [5], we present in this Appendix a somehow different 
procedure for justifying the reduction scheme presented in Section 7. 
Let the Taylor expansion of the 5— function: 

S(r - r') = E —J— X 'h ■ ■ ■ x 'in d h ■ ■ ■ x 'iJ( r )- 



n>0 



nl 



Applied to the current and charge distributions, it gives the following results: 



j(r,t)= j(r',t)5(r-r')d 3 x' = Yl 



■l) n 



n>0 



111 



X 



, h ...x' i J{r\t)d 6 x'd il ...d i J{v) 



P(r, t) = £ i-/- J < . . . x' in p(r', t)d 3 x'd h . . . d in 8{r). 

71>0 

The electric multipole moments are introduced by 



ii. ..i v 



x 



' n ...x' ln p(r',t)d 3 x' 



such that: 

P(r,t) = E ■ ■ ■ 9iJ(r) = E ^P (n) ||V^(r). 



nl 



nl 



n>0 n>0 

Considering the continuity equation for electric charge, we may write the equation 



ji = V'«?) + x 



,dp_ 
df 



Let 

(n) 

a) 



'n--- x 'ijidh---dij{r) 
x' ti ...x' l V'(x' l j)d 3 x'd ll ...d l J(r) + 
-Jx' i j-V'(x' il ...x'Jd il ...d i J(r) + 
-njx' h ... x'^x'ji^x'di, . . . d in 5{r) + 
-n x' h ... x' in _ x (x' i j ln - x' in ji) d 3 x'd h . . . d in 5{r) 



P (n+1) (t)||V n 
P (n+1) (t)||V" 



P (n+1) (t)||V" 



5(r) 
S(r) 
S(r) 



(n) , 

a) + 



P (n+1) (t)||V n 



S(r). 
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One obtains: 



» 



+ 



n 



n + 



71 + 1 



Y £ Hnk J < • • • {r x j) k d 3 x'd h . . . d in 5(r) 
S(r). 



P (n+1) (t)||V" 



The magnetic multipole tensor is introduced by: 



M, 



n 



, , X: . . . d ,(r'y. j) 4 dV 



and one may write 



» 



-e ak diM h „ An ^ lk d h . . . 9 in _ 1 5(r) + 



n + 1 



P (n+1) (t)||V 



5(r) 



with the result 

n>l 

+ E 

ra>l 



\n— 1 



^[M^HV"- 1 ! 5(r) 



*(r). 



Some invariance properties of the field are established in the following 

,(n) 



By introducing in the expression of j the tensor M/?l defined by the equation 



and using the relation £uk£i x kq - 



Ji = 



(M (n 



<)) + 



ji (M (n) 
(-l) n (n - 1) 



Mg ) 



-l) n (Yi - 1) 



rain 



n\n 

'Af[L (n) ]\\V n - 2 ].A5(r). 



hijiq + one obtains 

E ^^XVi J"- (n) ]^^i • • • d in ^S(r) 
X[^ n) ]\\V n ~ 2 } didi5(r) 



Considering the corresponding results for the potentials, the presence of a term like 
f(t)AS(r) in the expansion of j leads to a term in A or in $ which contains the 
expression 

A f(t-R/c) _ 

r c 2 r ' ' 

Therefore, regarding the contribution to the vector potential, the j expansion is 
equivalent to: 

(-l) n (n- I) 



j;=ji(M w -» 

;-i)>- 1; 



+ 



run 



m!2 

^ ;/ nine 



Af[l {n) ]\\V n - 2 ,5{r 



The last term in the precedent expression produces in A an additional term like V\l/(r, £) 
which, as one sees bellow, is a contribution from a gauge transformation of the potentials. 
Let the transformation 



P; 



«l...i n _l 



" 1 ^,., n _ 1 [LW] 



c 2 n 2 
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By substituting Pj 1 ...i n _ 1 in the j and p expansions by 

n — 1 • 



«i...i n _i 



p' . _|_ — 2|\| 

r n...t n _ 1 ^ c 2 n 2 in-i- 



one gets: 



p(n-l) 



(n) 
p/(n-l) 



-l) n (n - I) 



n\n 



' N[\S n) ]\\V n - 1 } didi8{r) 



(A.1) 



-l) n fn- 1' 



nine' 



"^'[L (n) ]||V n - 1 l <5(r). 



It is easy to see that the potentials ($, A) corresponding to the new densities 
(P, j) with 



) _> p'(n-l) 



differ from the original one by a gauge transformation. 

We point out that the physical equivalence of the two distributions (p, j) and 
(p, j) is true only for the field in the exterior of the domain T> including the support of 
the charges and currents. 

• Let 

J 71—1 

P(L)ii...j„ — Pti...i„ — ~~ £ i i i„q'^'a 1 „,i„_ 1 (][^ 



A=l 



It is a simple matter to see that the transformation — > Pff does not alter the 
density p but 



Ji = J* 



'p( n ) 



nln 



-euA[M il ..^AL}\\V n - 2 \ q 5(r) 



The difference between ji[P^] and ji^P^ — > P[n] m ay be eliminated by the 



transformation 



and so 



Ti — 1 • 

M,,.,.. . — > M' ll ..., n _ 1 M il ... in _ 1 + — ^..^[L] 



Let the transformation (J37j) of the magnetic tensor 



M ( 5)u... 4 „ = Mii...i„ - E ^ 2 A* 3 ...US (n) ] 

£»(<) 



where S*-"- 1 is a fully symmetric tensor and the operator A is used as in equation (J35)l for 
detracing an arbytrary fully symmetric tensor. 

Introducing the tensor in the expansion of j we may write 



n-1 



j, = ji(M'-' — M( 5 ' } j + - — - — e Uk di ^ 1 i 2 A i3 ...i n _ lfc [S (n) ]9 il . ..d in ^S(r) 

D(i) 
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A being a fully symmetric tensor. In the sum there are n — 1 terms with 8^, j = 
1 . . . n— 1 but Ski j di j = dk and e^didk = such that these terms give a null result. In the 
remaining terms we have the factors like 8^, j,q — 1 . . . n — 1 resulting 5i j i q di j di q = A. 
The number of these terms is C\_ x = (n — l)(n — 2)/2, such that 



/ M (n) , M ("K | 1 (n-l){n-2) (n) 



,M(r) 



iM n) - m$) + ( ~ irl(n 2 ~ 1)(n ~ 2) {V x [A[SW]||V n - 3 ]}. A5(r 



This relation allows the introduction of an equivalent current density regarding the 
vector potential: 

(-l) n - 1 (n-l)(n-2) 



ji" =ji{M 



(n) 



( n ) \ i 
M (S ) + 



{V x [A[S (n) ]||V n - 3 ]}.5(r). 



{S)j 1 2n!c 2 
The suplimentary term is eliminated by the transformation 

M( n - 2 ) -> M ,(n ~ 2) = M (n ~ 2) + ^~ 2 A[S (n) ] 



2nc 2 



such that 



M (n-2) 



M 



M 



(n) 
(5) 
/(n-2) 



Let 



r>(») 

One gets 



ri! 

vn— 1 / 



j, = j,(P ( "»^P? , ) + 
= ji (p.^pW) + (-^-'("-i) 



r>(0 



A[S^]||V n - 2 l^(r) 



-l)(n- 2) 
2n! 



A[S (n) ]||V n - 3 l.A5(r)' 



So, we obtain an equivalent current density 



)(»») 



+ 



Ji — Ji {' ~^ 1 5 

(-l)™- 1 ^- l)(n-2) r 
+ 2d^2 
For the charge density we have 



n! 



A[S (n) ]||V n - 2 9i<J(r) 



A [S (n) ]||V™- 3 l.(5(r). 



p = p (P (n) - pg°) + ^r 1 E ^a^.us^]^ . . . ^*(r). 

n - D(i) 

Because in the last sum there are C 2 terms like 8{a with equal contributions, one obtains 



P = P (P W 



+ 



;-l) n n(n- 1) 
2d 



A[S 



(n)i I 



A5(r) 
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P = P (P (n) - pP) + 



2n!c 2 



A[sW]||v n " 2 5(r) 



Let us define 



p"(n-2) = p (n-2) « -A[S( n )l 

2nc 2 1 J 



It is easy to see that 



I p (n-2) p "(n-2) 



n! 



and 



,(n) __, P W 



P = P 



p (n-2) _^ p "(n-2) 

So, the densities (p',j') and 

f / p (n) ^ p(«) 



,(n) _, p (") 



P I p(n-2) _^ p "(n-2) I > -M p (n-2) _^ p "(n-2) 



are equivalent because the corresponding potentials differ only by a gauge 
transformation. 

Appendix B. The Reduction Scheme 



(jjl, e) = (l, 2) 



p(2) 



7?(2) 



Note "P = P - -P 8 - 
no extra- gauge changes for A and $ for P 1 - 2 - 1 



p(2) 
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(jjl,6) = (2,3) 



M< 2 > 



N (2,l) 









p(3) 




n 


(i) 





p(3) 



Eq.(I),n=2 



Eq.(IV),n=3 



Note:T, = IN? ' 1} - |fl, 



P (i) _ ^Af[M^] = P« 



PW - £AT[M< 2 >] 



4c 2 " 



(2,1) 



6<r 



(1) 



P (1) - it 



(1) 
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(/*,£) = (3,4) 



M< 3 > 



N (3,l) 



M 



(3) 

(sym) 



A 



(i) 



p(4) 



n 




N 



(3,2) 



P (2) - ^M ym [M( 3 )] + ^ApW 



Note. N ( i:i, } = N (3,1) = traceless 



sym 



| N (3 ' 1} - _ i n , 



Transf. Il( 2 ) -> ft 2 ) 

produces only a 
gauge transformation 



Eq.(I),n=3 



Eq.(III),n=3 



Eq.(IV),n=4 



Eq.(II),n=3 



A[n( 2 )] 




p(2) 






Mp (4) ] 


= ^(2) _ _^t (2) 




M (1) + ^A[Mg)j 
M« + ^A« 




M (1) + £A[Mg> 
i 

18c 2 



6c 2 ' 

1 A/" 2 fM (3) l 



syml 



M« + ^A« 

T^^ (3 ' 2) 
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Gu,e) = (4,5) 



N (4,l) 



A (2) 



p(5) 



n 



(3) 





p(5) 




T>(3) _ 3 f^ 4 ' 1 ) , 3 jj(3) 
' 16c 2 1X1 ^ 10c 2 11 



N (4 ' 2) 



T>(3) _ ^|\| (4,1) + _3_f[(3) 
' 16c 2 'Vm ^ 10c 2 11 



a[n<$], A[n(^ 



~(4,1) ~(3) 

t>(3) 3 N v = p( 3 )- 

' 16c 2 1X1 ^ 10c 2 11 ' 



^T (3) 

c 2 



A^( 2 ) + £A(2) - 



(4,2) 



N 



(4,3) 



A// (2) _j_ 1 A (2) _ 1 |\| {4,2) 
JVI ^ 4c 2 YV 24c 2 ™sym 

A... 



Eq.(I),n=4 




Eq.{III),n=A 



Eq.(IV),n=5 



Eq.(II),n=4 



Eq.(IV),n=3 



Eq.(I),n=2 



no effect 



= M<® + 





^A[P {5) ] - P (1) - 

3^A[NiS + iA[n(3)] 



l_ t (l) _ 

„2 1 



p(1) -^ t(1) -3 2 ^ A ( A/ '^[ M(4) ]] + 

^A[A[P( 5 )]] + ^>W] = 
P (1) - £t (1) - ^A[NgS] + 



i 

20C 



-ah": • w^n 



.-,(4,3) 



□ 
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Appendix C. Results 



Table I 





(ftO = (1,2) 


(//,£) = (2, 3) 


(//,£) = (3, 4) 




P<« : (p) 


pd) _ 4t (1) 

C 


pd) _ 4t (1) 

c 


p<» 


17(2) . 17 p lp r 

' • ' »j r u 3 r qq u v 


p(2) 


p(2) _ JL_t (2) 


p<») 




77(3) 


77(3) 


pW 






77(4) 


pO» 








M (1) 


M (1) : (m) 


M« 


M« + 1 A(i) - T ----N (3 ' 2) 


M ,2) 








M M 








M ,4 > 









Table II 





(/<,£) = (4, 5) 


pW 


p (1 » - - ^ A [NSSl + ^A[nM] + ^ n* 4 ' 3 ' 


pM 


77(2) _ J_j {2) 


pW 


p(3) - ^t (3) 


pW 


77(4) 


p(=> 


77(5) 


M W 




M ,2) 


„ . ~(2) ~(4,2) 

-M (2) + M -5^N 


m' 3 ' 




m' 4 ' 





In the above formulas it was introduced the electric toroid moment of undefined 
order, as it results from the algorithm presented before: 

Tin) = n .(n+1,1) _ n ~M 
(n+1) 2 2(n + 2) 
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